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ABSTRACT: An analytical and experimental investigation was carried out on G40- 

800/R6376 graphite epoxy laminates to evaluate the Edge Crack Torsion (ECT) test as a 
candidate for a standard Mode HI interlaminar fracture toughness test for laminated composites. 
The ECT test consists of a [90/(±45)3/(+45) 3 /90] s laminate with a delamination introduced by 
a non-adhesive film at the mid-plane along one edge and loaded in a special fixture to create 
torsion along the length of the laminate. Dye penetrate enhanced X-radiograph of failed 
specimens revealed that the delamination initiated at the middle of the specimen length and 
propagated in a self similar manner along the laminate mid-plane. A three-dimensional finite 
element analysis was performed that indicated that a pure Mode HI delamination exists at the 
middle of specimen length away from both ends. At the ends near the loading point a small 
Mode II component exists. However, the magnitude of this Mode II strain energy release rate 
at the loading point is small compared to the magnitude of Mode HI component in the mid- 
section of the specimen. Hence, the ECT test yielded the desired Mode HI delamination. The 
Mode in fracture toughness was obtained from a compliance calibration method and was in 
good agreement with the finite element results. Mode II End-Notched Flexure (ENF) tests and 
Mode I Double Cantilever Beam (DCB) tests were also performed for the same composite 
material. The Mode I fracture toughness was much smaller than both the Mode n and Mode in 
fracture toughness. The Mode H fracture toughness was found to be 75% of the Mode in 
fracture toughness. 

KEY WORDS: composite materials, fracture toughness, delamination, compliance 
calibration, three-dimensional finite element analysis, strain energy release rate 


Nomenclature 

a delamination length for DCB and ENF specimens, m 
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a delamination length for ECT specimens, m 
Aa virtual crack length, m 

b specimen width for DCB and ENF specimens, m 
B specimen width for ECT specimens, m 
C loading point compliance, m-N' 1 

5 loading point displacement, m 

5xi relative crack face displacement in x-direction at node i and i', m 

5 y i relative crack face displacement in y-direction at node i and i', m 

8zi relative crack face displacement in z-direction at node i and i', m 

Eh Young's moduli, i=l, 2, 3, GPa 

Fxj nodal force in x-direction at node j 

Fyj nodal force in y-direction at node j 

F Z j nodal force in z-direction at node j 

Gy shear moduli, ij=12, 13, 23, GPa 

G strain energy release rate, J-m' 2 

Gi Mode I strain energy release rate component, J-m' 2 
Gn Mode II strain energy release rate component, J-m' 2 
Gm Mode HI strain energy release rate component, J-m' 2 
Gjc Mode I fracture toughness, J-m' 2 

Gnc Mode II fracture toughness, J-m' 2 

Gjnc Mode III fracture toughness, kJ-m' 2 

h specimen thickness, m 

H ply thickness, m 

L effective ECT specimen length, m 

L ' specimen length, m 

21 span of ENF specimen, m 
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m linear regression parameter for ECT tests 
m ' linear regression parameter for ENF tests, N'^m' 2 
n linear regression parameter for DCB tests 

Vij Poisson’s ratio, ij=12, 13 and 23 

P transverse load, N 

P c transverse load at delamination growth, N 
W moment arm length for ECT specimens, m 

x, y, z Cartesian coordinates 

Introduction 

Recently the Edge Crack Torsion (ECT) test was proposed as a candidate standard 
Mode HI fracture toughness test for laminated composites [1], The ECT test consists of a 
[90/(±45) 3 /(+45) 3 /90 ] s laminate with a delamination introduced by a non-adhesive film 
inserted at the mid-plane along one edge during manufacture and loaded in a special fixture to 
create torsion along the length of the laminate. Previous analytical investigations of the ECT 
test [2,3] indicated that the failure mode is Mode HI in the portion of the specimen away from 
both ends where the load nose and supports contact the specimen. However, experimental 
evidence of Mode HI delamination growth from the insert is needed to justify the validity of the 
ECT test. In addition, to understand the boundary effects, a full three-dimensional finite 
element analysis is required. 

The first objective of this paper was to understand and document the nature of 
delamination growth in the ECT specimen. The second objective was to develop a three- 
dimensional finite element model of the ECT specimen to accurately evaluate the individual 
mode and total strain energy release rate along the delamination front, taking into account the 
three-dimensional effects neglected by previous analytical investigation [2,3]. To this end. 
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G40-800/R6376 graphite/epoxy laminates were manufactured by Ciba Composites and were 
tested monotonically to failure in the ECT fixture. A three-dimensional finite element analysis 
of a typical ECT specimen was also performed. To compare the magnitude of Mode in fracture 
toughness to the Mode H fracture toughness. Mode H tests were performed using End-Notched 
Hexture (ENF) specimens. In addition. Mode I Double Cantilever Beam (DCB) tests were also 
performed for a complete comparison. In the paper, the test results for the ECT specimens, the 
ENF specimens and the DCB specimens were presented first, followed by the three- 
dimensional finite element analysis of a typical specimen. 

Experimental Investigation 

Specimen Configuration 

The ECT laminate specimen shown in Figure 1 was an 89 mm long (L') by 38 mm 
wide (B) (3.5 in. x 1.5 in.) rectangular plate, where the span (L) and the moment arm (W) 
were 76 mm (3 in.) and 32 mm (1.25 in.), respectively. The average specimen thickness (h) 
was 4.9 mm (0.193 in.). The ECT specimens were cut from a [90/(±45)3/(+45) 3 /90] s 
laminated panel with embedded 13 |im non-adhesive Polytetrafluoroethylene (PTFE) films 
embedded at the mid-plane along one edge to introduce an edge delamination crack of length, a. 
The panel was made of G40-800/R6376 graphite/epoxy tape plies. The fiber aerial weight 
(FAW) and fiber density (FD) were 190 g nr 2 and 1.80 g em 3 , respectively. The fiber volume 
fraction (V f ) for the ECT specimens was calculated to be 60 %, (percent V f =(Number of 
plies)*(FAW)/[h (FD)]100). C-scans were used to determine the insert location in the panel 
which helped when cutting the panel into specimens. Due to the unstable delamination growth 
of the ECT specimen [1], a series of specimens with initial delamination lengths varying 
between 0 to 19 mm (0.75 in.) were made to generate a compliance calibration curve. Three 
groups of specimens, with six replicates each, were cut from a single panel. Within each 
group, the insert lengths (initial delamination lengths, a) were 0, 2.5, 5, 11, 16 and 19 mm, 
respectively. Note the delamination length a was measured from the free edge to the end of 
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insert while the delamination length, a, was measured from the loadline to the end of insert for 
both DCB and ENF specimens. 

Experimental Procedure 

The ECT specimen was placed in the ECT test fixture shown in Figure 2. The fixture 
provided supports at three points, and the transverse load was introduced at the fourth point by 
compression of the loading pin. The loading pin was loaded in the test machine with a constant 
stroke rate of 1.5 mmmin-i. During the test, the specimen was subjected to end torsion 
induced by a pair of counteracting moments at a distance L=76 mm (3 in.) apart Each moment 
was generated by a pair of forces with a moment arm W=32 mm (1.25 in.). The applied load P 
versus load-point displacement 8 was recorded on an X-Y plotter for each specimen. 

Test Results 

A typical load-displacement (P-8) curve was shown in Figure 3. The curve was nearly 
linear up to the point A. After this point, the delamination grew from the insert and the stiffness 
decreased until a load drop was observed at P c . A line tangent to the initial linear portion of the 
loading curve was drawn to determine the initial compliance C=8/P. Initial deviation from 
linearity due to seating of the loading fixture was ignored. The critical load for delamination 
propagation was defined as P c . The critical load P c along with the initial delamination length 
was used to calculate the critical total strain energy release rate G c . 

Dye penetrate enhanced X-radiographs of the ECT specimens, taken after the load drop 
occurred, were used to identify the delamination growth from the insert. Figures 4 (a) and (b) 
showed two different shapes of the delamination front for specimens with insert widths of 1 1 
mm and 2.5 mm, respectively. In Figure 4 (a), the delamination initiated at the end of the insert 
in the middle of the effective specimen length and propagated in a self similar manner along the 
laminate mid-plane between the two 90° plies. Hence, the desired Mode III delamination 
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growth pattern was achieved in this specimen. 


However, for the specimen with the 2.5 mm insert, shown in Figure 4 (b), the 
delamination grew off center of the effective specimen length. An examination of the edge of 
the specimen indicated that delaminations grew at the 90/45 and the 45/-45 interfaces. No 
delamination growth was evident between the two mid-plane 90° plies where the insert was 
placed. Thus Mode III type delamination growth was not achieved in this specimen. This 
behavior could be attributed to the fact that the insert length was too small. As observed in 
Ref. 2, the strain energy release rate dropped sharply from a plateau value to zero as the 
normalized initial delamination length (a/B) reduced from about 20% to zero under the same 
twist. In Figure 4(b) the normalized initial delamination length was (a/B)=0.07. This value 
appears to be too small to have enough delamination driving force at the insert front to initiate 
Mode in type delamination growth. Consequently, specimens with insert length of 2.5 and 5 
mm were discarded. 


Mode m-Compliance Calibration Method 

When the critical strain energy release rate is established experimentally from measuring 
the change in compliance of the cracked body as a function of crack length, the data reduction 
method is usually referred to as the compliance calibration method. As shown in Ref. 4, the 
strain energy release rate can be calculated from the change of compliance with crack length of 
the delaminated specimen as 


G 


p^ac 

2L 3a 


where 


3C 

da 


is the rate of change in compliance as the crack grows. 


( 1 ) 
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Two procedures are often adopted in the experimental evaluation to approximate the 
derivative of C with respective to a. In the first procedure, the delamination is extended by a 
small increment Aa and the change in compliance AC is measured. Then the ratio (AC/Aa) is 
used to approximate the derivative in Equation 1. In the second procedure, the compliance C is 
measured by varying the initial crack length a by moving the specimen in the test fixture and 
applying a small load to establish the specimen compliance without growing the crack. A plot 
of C versus a is then generated. A least squares curve fit is applied to determine C as a function 
of crack length. The derivative of compliance as a function of a can then be substituted into 
Equation 1. 

Unfortunately, for the ECT configuration it is not possible to generate a compliance 
curve from a single ECT specimen using either method because of (1) the unstable 
delamination growth and (2) it is not possible to shift the specimen in the fixture and change the 
initial delamination length, a. Hence, a number of specimens with different insert lengths were 
tested to obtained a plot of compliance as a function of delamination length. First, the 1/C (the 
inverse of compliance) values were plotted as a function of normalized delamination length a/B 
as shown in Figure 5. Then, a least square straight line fit to the data was determined. The 
linear relationship between (1/C) and (a/B) represented by the straight line can be expressed as 



where m and A are constants. The Mode III fracture toughness, G mc , can be obtained by 
differentiating Equation (2) with respect to delamination length a and substituting the derivative 
of C with respect to a into Equation (1), yielding 
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(3) 


Gmc ~ 


mPfe 


2LB 

'l-mf-YI 




For the data presented in Figure 5, m and A were obtained as 0.9 and 415 kN-nr 1 . The 
Mode ID fracture toughness G mc of the individual specimens can be calculated using measured 
values of C and P c from the individual specimens and the value of m from the curve fit 
described above. The critical load P c , the inverse of the compliance 1/C, and the fracture 
toughness were presented in Table 1 along with the specimen thickness and initial delamination 
length. A more accurate delamination length measurement was made after the test by 
completely opening the delamination and observing the delamination surfaces. The end of the 
insert could be readily distinguished from the propagated delamination on the delaminated 
surfaces. For specimens where the insert front was not straight, the delamination length was 
determined from the average of the insert length measurements at 6 or more equally spaced 
locations along the insert front. If these measurements showed large variations, the results 
should be discarded. Specimens with insert lengths measured as 11 mm, 16 mm and 19 mm 
before the test were opened after the test. The insert front was fairly straight in these specimens 
(within 1 mm) and hence no results were discarded. The measured insert length after the test 
was within 1 mm of the initial measurement for each of these specimens. For the nine G40- 
800/R6376 graphite epoxy specimens tested, the mean Mode HI fracture toughness, G mc , was 
1419 J nr 2 with a coefficient of variation of 8.2%. 

Comparison with Other Delamination Fracture Modes 

Mode II Fracture Toughness 

In order to compare the Mode HI fracture toughness to the Mode II fracture toughness 
of the G40-800/R6367 graphite epoxy material, End-Notched Flexure (ENF) Mode II tests 
[5,6] were also performed. The ENF specimen configuration, shown in Figure 6 (a), consisted 
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of a unidirectional laminate with a mid-plane embedded insert at one end subjected to 3-point 
bending loading. The total specimen length was about 250 mm, width b was 25 mm and the 
thickness h was 4.3 mm (24 unidirectional plies). The length between the two supports ( 21 ) is 
100 mm. A 13 pm non-adhesive PTFE film was placed at the laminate mid-plane during 
manufacture. The insert length was 60 mm for each specimen. The ENF apparatus is shown in 
Figure 6 (b) where a spring-loaded displacement gauge is mounted under the specimen as 
shown. 


The specimens were loaded and unloaded in the elastic range first to obtain compliance 
with respect to initial delamination lengths of a=0, 15, 20, 25, 30, 35, and 40 mm. This was 
achieved by keeping the supports fixed and sliding the specimen in the rig to achieve different 
initial delamination lengths. Note that sliding the specimen produced different distances 
between the right roller support and the tip of the insert (a in Figure 6(a)). Specimens were 
loaded in displacement control loading with a constant stroke rate of 0.5 mm/min. After the 
compliance measurements, a least squares linear regression of the following form was 
performed. 


C = C 0 + m' a 3 


(4) 


The results of the compliance calibration tests for each specimen were presented in Table 2 
along with the corresponding value of m’. 


The specimen was then positioned in the test fixture so that aAf=0.5, and loaded until the 
delamination started to grow from the insert. The Mode II fracture toughness was evaluated 
based on the compliance calibration [6] as 


»Hc 


_ 3m' a 2 Pg 

2b 


(5) 
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The maximum load and the initial delamination length was used in Equation 5 to 
calculate G nc . A total of five specimens were tested. The test results were presented in Table 3. 
One specimen was completely separated after the test to verify the insert length. The mean 
Mode II fracture toughness G nc was 1061 Jnr 2 with a coefficient of variation of 9.8%. This 
value of Mode II fracture toughness was 75% of the Mode III fracture toughness for this 
material. 

Mode I Fracture Toughness 

For a complete comparison. Mode I Double Cantilever Beam (DCB) tests were also 
performed in a manner largely following the ASTM Test Method for Mode I Interlaminar 
Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites (D 5528- 
94a). A 13 (im thick non- adhesive PTFE film was placed in the mid-plane of the laminate 
during fabrication to introduce initial delamination to the specimens that were subsequently cut 
from the laminate. The DCB test specimens were 4.3 mm thick (h) (24 ply unidirectional 
laminate), 12.5 mm wide (b) and 190 mm long (L'). The 0° fibers were oriented along the 
length of the specimens. A schematic of the DCB specimen was shown in Figure 7. For each 
specimen, the insert length was about 45 mm from the specimen end or 38 mm from the load 
line. To introduce load to the specimen, loading blocks of 13 mm x 13 mm x 13 mm were 
bonded to the specimen end using an adhesive. A total of five specimens were tested. 

The DCB specimens were loaded in a test machine at a displacement rate of 0.5 
mm-min 1 to open and propagate the cracks. The specimen edges were painted with a white 
typewriter correction fluid to help identify the delamination front during the tests. Initial 
delamination growth from the insert for each specimen was determined from the sudden load 
drop of the initial linear load-displacement trace. As the delamination grew further, the location 
of the crack front was marked on the specimen edge at a crack growth interval of about 6 mm. 
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For each marked crack front, the corresponding location on the load-displacement trace was 
also recorded for later data reduction. 


To calculate the Mode I fracture toughness, Gic, the compliance calibration method was 
used for data reduction [7]. A least squares fit of the log-log plot of specimen compliance C= 
8/P as a function of delamination length was performed to obtain a slop n where C=K-a n 
(where K is a constant). The Mode I fracture toughness Gic was determined from. 


Gic - 


nPS 

2ba 


( 6 ) 


The initial Gic values were calculated from n along with the P, 5 and a values 
associated with the delamination onset from the insert. These Gic values and the n value 
obtained were given in Table 4. The average initial Gic value was 322 J-nr 2 with a coefficient 
of variation of 7.6%. 


The comparison of various modes of fracture toughness from the DCB, ENF and ECT 
tests was shown in Figure 8. The DCB Mode I fracture toughness was about 30% of the ENF 
Mode II fracture toughness and 23% of the ECT Mode m fracture toughness for the G40- 
80Q/R6376 graphite epoxy studied. 


Analytical Investigation 

Finite Element Model 

A finite element model for ECT specimen C5 in Table 1 with dimensions 89x38x4.8 
mm (3.5x 1.5 x0.189 in.) was developed using 34,048 eight-node solid elements (SOLID45). 
The ANSYS finite element analysis software package (Version 5.1) was used for the analysis. 
A finite element mesh of the ECT specimen is shown in Figure 9 along with a cross section 
view of the detailed finite element mesh at the delamination front and an enlarged view of the 
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mesh at a specimen comer. Each ply thickness (H) was assumed to be 0.171 mm (0.00675 
in.). The initial crack length of 19 mm (0.75 in.) was assumed. The crack plane was modeled 
with double nodes of zero distance connected by 4669 frictionless point-to-surface contact 
elements (CONTAC49) to prevent interpenetration of the two delamination surfaces. The 
specimen length was modeled with 28 equally spaced elements. The specimen width was 
modeled with 38 elements distributed symmetrically about the half width. The element size at 
the half width where the crack tip was located was 0.171 mm (0.00675 in.), and increased 
linearly to 3.175 mm (0.125 in.) at the two edges. The thickness of each ply was modeled by 
one element except for the two 90° plies on either side of the mid-plane. These plies were each 
modeled with three equal sized elements through-the-thickness. Therefore the mesh size of the 
elements surrounding the delamination crack tip was l/3HxH (see Figure 9b). 

The lamina elastic properties were obtained from tension tests of unidirectional and 
±45° laminate coupons and were presented in Table 5. In the analysis, G 23 = G 12 and v 23 = 
V 12 were assumed because G 23 or v 23 are much more difficult to determine experimentally. 
According to the analyses presented in reference 2, these assumptions result in conservative 
estimation of the critical strain energy release rate compared with other G 23 values 
(0<G 23 <G 12 ). The effect of these assumptions on the prediction of the strain energy release 
rate components were studied by the present three-dimensional finite element model and the 
results were presented in the Appendix. As shown in the Appendix the assumptions of G 23 = 
G 12 and v 23 = v 12 resulted in slightly lower Mode HI strain energy release rate component 
than the case where G 23 =V f -G 12 and v 23 =(E 22 /2G 23 -l) were used. For G 23 values differing 
by 40%, the differences in the average total strain energy release rate and Mode in component 
were 9% and 10%. The effect of G 23 and v 23 on the strain energy release rate values may be 
considered small. The three-dimensional finite element analysis confirmed the two-dimensional 
analytical predictions on the G 23 effect on strain energy release rate. 
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The model was supported in z-direction at three nodes where the nodal locations were 
equivalent to the position of the three supporting pins on the ECT test fixture. The remaining 
two translational rigid body motions were also prevented by constraining the degree of 
freedoms along x-direction and y-direction at two of these supporting nodes. A critical test load 
of 1045 N (235 pounds) for specimen C5 (corresponding to a load point displacement of 4.63 
mm) was applied at a node that was located 3.175 mm (0.125 inch) along the specimen length, 
and 3.175 mm (0.125 inch) in from the edge, from the specimen comer. Because the ECT 
model contained contact elements, the analysis required a geometrically nonlinear solution to 
update the deformed geometry during the solution process. However, contact was only found 
in the local vicinity near the loading points as shown in Figure 10, and hence, did not influence 
the calculated strain energy release rate components. 

Virtual Crack Closure Techniq ue 

The virtual-crack-closure-technique (VCCT) [7] was used to calculate the strain energy 
release rate components for Mode I, Mode II and Mode IE. The VCCT uses the stresses/nodal 
forces ahead of the delamination front and the displacements behind the front to determine the 
G components. A simple formula for calculating the strain energy release rate components for 
two-dimensional four- node elements was given by Rybicki and Kanninen [8]. This formula 
uses the forces at the crack tip and the relative displacements of the crack faces behind the crack 
tip to calculate the G components. The strain energy release rate components can be calculated 
from the work required to close the delamination by one element length Aa (see Figure 1 1 (a)) 
as 


Mg 

2Aa 


Gr = -S-B. 


Vy 

2Aa 


r _ yj yi 

u n - 


(7) 

( 8 ) 
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where F y j and F z j are the forces in the y-direction and z-direction at node j, 8^ and 8 zi are the 
relative crack face displacements between nodes i and i', located at a distance Aa behind the 
crack tip in the y-direction and z-direction, respectively. 


Extension of these formulas to three-dimensional eight-node solid elements is straight 
forward. The strain energy release rate components are calculated (see Figure 1 1 (b)) from 


Fzi^zi 

G x =— 

2AaAx 


g IsAl 
u 2AaAx 


r* _ xjxi 
'Jffi 


F-A 

2AaAx 


(9) 


(10) 


( 11 ) 


where Ax is the sum of one-half the element lengths on either side of node j in the x-direction. 

Nymepcal Resets 

The typical distribution of strain energy release rate components along the length of the 
delamination front was shown in Figure 12. The Mode I component was negligibly small along 
the entire delamination front. A finite Mode II component was observed near the loading 
points, but diminished away from the loading point The Mode HI component raised from zero 
at both ends and reached a plateau about 20% into the specimen length from either support The 
plateau value of the Mode in component was much higher (by a factor nearly equal to 3) than 
the maximum Mode II value observed at the loading point. Note that the Mode n fracture 
toughness was about 75% of the fracture toughness from the ECT tests in previous section, 
hence, the magnitude of the Mode H component at the loading point will not reach the Mode H 
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fracture toughness when the magnitude of Mode HI component in the mid-section of the 
specimen reaches the Mode IE fracture toughness. Hence, a Mode El delamination will initiate 
in the ECT specimen. 

The three-dimensional finite element analysis shows that the delamination will grow in 
pure Mode El at the middle of the specimen length away from both ends. This is in agreement 
with the experimental observation from the dye penetrate enhanced X-radiograph shown in 
Figure 4 (a). An average strain energy release rate component can be calculated by integrating 
that component over the effective region along the delamination front and then divided by the 
effective length of the ECT specimen. The average total strain energy release rate and the Mode 
TTT component computed from the finite element model were 1349 J-nr 2 and 1239 J-nr 2 , 
respectively, which indicated the critical total strain energy release rate integrated along the 
delamination front was 92% Mode IE. The critical total strain energy release rate for the ECT 
specimen can be used to approximate the Mode IE fracture toughness. The finite element result 
for the average total strain energy release rate agreed very well with the mean test result of 1419 
J-nr 2 determined using compliance calibration. 

Concluding Remarks 

An analytical and experimental investigation was carried out on G40-800/R6376 
graphite epoxy laminates to evaluate the Edge Crack Torsion (ECT) test as a candidate standard 
Mode IE interlaminar fracture toughness test for laminated composites. Dye penetrate enhanced 
X-radiographs of failed specimens revealed that the delamination initiated at the middle of the 
specimen length and propagated in a self similar manner in the laminate mid-plane. A three- 
dimensional finite element analysis was performed that indicated that a pure Mode El 
delamination will exist at the middle of specimen length away from both ends. A Mode E 
component was also shown to exist near the ends where the loading points are located. 
However, the magnitude of the Mode E strain energy release rate at the loading point was 
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much smaller than the magnitude of Mode in component in the mid-section of the specimen. 
Hence, the failure mode was the anticipated Mode D3 delamination. The total strain energy 
release rate calculated from the finite element model, with the applied load (P) equaled the 
critical load (P c ) observed in the experiment, was in good agreement with the fracture 
toughness evaluated from the compliance calibration method for the ECT specimen. The mean 
Mode IE fracture toughness from nine specimens was 1419 Jnr 2 with a coefficient of 
variation of 8.2%. Mode I DCB and Mode n ENF tests were also performed on the same 
composite material. The mean Mode D fracture toughness from five specimens was found to be 
75% of the ECT Mode HI fracture toughness. The average Mode I fracture toughness from five 
DCB specimens was about 30% of the ENF Mode H fracture toughness and 23% of the ECT 
Mode HI fracture toughness. 

Acknowledgment 

This work was performed while the first author was a National Research Council Research 
Associate with the Army Research Laboratory Vehicle Structures Directorate at NASA Langley 
Research Center. 

Appendix 

Shear Modulus G 2 3 Effect on the Calculation of Strain Energy Release Rate 

The purpose of this appendix was to study the effect of G 2 3 and v 2 3 on the strain 
energy release rate components for the ECT specimen under consideration. The three- 
dimensional finite element model for specimen C5 shown in Figure 9 was re-analyzed under 
the same loading conditions with the material properties G 23 =V f -G 12 =2.91 GPa and 
v 23 =(E 22 /2G 23 -l)=0.54 while retaining the other material properties in Table 5. The 
distributions of the strain energy release rate components along the delamination front were 
plotted in Figure 13 along with the results shown in Figure 12 for comparisons. Case 1 
represents the results from the three-dimensional finite element analysis by using the 
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assumptions of G 23 =G 12 =5.15 GPa and v 23 =v 12 =0.325 while case 2 represents the result 
with the properties G 23 =2.91 GPa and v 23 =0.54. Figure 13 shows that the distributions of the 
individual mode strain energy release rate are similar for these two cases. The Mode III 
component was higher for case 2 than case 1 by about 10%. The Mode II component was 
nearly the same throughout the effective length of the ECT specimen in both cases except at the 
loading and support locations where it was lower in case 2 than in case 1. The average total 
strain energy release rate and Mode in component for both cases are summarized in the 



Average Strain Energy Release Rate J-nr 2 


Total 

Mode HI 

Case 1 

1349 

1239 

Case 2 

1480 

1371 


The average total strain energy release rate and Mode HI strain energy release rate component 
obtained from the two cases differed by about nine and ten percent for G 23 values differing by 
40%. Hence the effect of G^ and v 23 on the strain energy release rate calculation appears to be 
small. The assumptions on G 23 and v 23 in case 1 give conservative predictions on the total 
strain energy release rate and the Mode HI component. These observations are in agreement 
with the two-dimensional analytical predictions given in reference 2. 
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Table 1. Test Results of the ECT Specimens 


Specimen No. 

Thickness (mm) 

a* (mm) 

1/C (kN-m-i) 

JPcOO 

Gmc (J-m-2) 

A3 

5.1 

11 

347 

1535 

1411 

A4 

5.0 

16 

267 

1165 

1253 

A5 

4.7 

19 

226 

1045 

1341 

B3 

4.9 

11 

316 

1432 

1352 

B4 

4.8 

16 

267 

1343 

1657 

B5 

4.8 

19 

211 

1050 

1448 

C3 

5.0 

11 

317 

1459 

1397 

C4 

4.9 

16 

270 

1236 

1389 

C5 

4.8 

19 

200 

1045 

1519 


♦Specimens with 2.5 and 5 mm delamination lengths were discarded. 


Table 2. Compliance Calibration Tests Results for ENF Specimens 


Compliance C for Specimen 1-5 (m-N _1 x 10' 6 ) 

a (mm) 

1 

2 

3 

4 

5 

0 

1.171 

1.219 

1.256 

1.142 

1.199 

15 

1.271 

1.285 

1.328 

1.179 

1.262 

20 

1.313 

1.342 

1.428 

1.291 

1.319 

25 

1.385 

1.442 

1.542 

1.371 

1.419 

30 

1.542 

1.613 

1.727 

1.536 

1.548 

35 

1.799 

1.856 

1.970 

1.767 

1.797 

40 

2.018 

2.150 

2.322 

2.048 

2.087 

Co 

1.198 

1.225 

1.277 

1.151 

1.204 

m' [N^-nr 2 ] 

131.1 

145.1 

163.6 

141.7 

137.1 
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Table 3. Test Results of the ENF Specimens 

Specimen No. 

Thickness (mm) 

8 C (mm) 

P C (N) 

Guc (J-m- 2 ) 

1 

4.4 

2.223 

1557 

1172 

2 

4.3 

2.051 

1379 

1018 

3 

4.3 

2.140 

1374 

1143 

4 

4.4 

2.000 

1423 

1061 

5 

4.3 

2.096 

1343 

912 


Table 4 Gic Results Obtained from DCB Tests Using Compliance Calibration Method 



(C=Ka n , 

n=2.744, Gic = n-P-8/(2-b a)) 


Specimen No. 

a. Crack Length (mm) P, Load (N) 8 , 

Displacement (mm) Gic (J-m* 2 ) 

1 

37.8 

43.9 

2.36 

296 

2 

37.1 

47.9 

2.44 

340 

3 

38.3 

42.1 

2.44 

289 

4 

37.4 

46.2 

2.51 

335 

5 

37.9 

51.7 

2.36 

349 


Table 5. Elastic pr operties of G40-800/R6376 graphite epoxy com posite laminate 

Properties 

E n = 151 GPa 

E 22 = E 33 =8.97 GPa 

G 12 = Gjj =5.15 GPa 

v i 2 = v i 3 =0-325 
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Figure 2. An ECT Specimen in the ECT fixture with load P and displacement 8 applied to the loading pin. 






Delamination Front 



Delamination Front 


Insert 



(b) Insert Length of 2.5 mm 


Figure 4. Dye penetrate enhanced X-radiographs of failed ECT specimens for insert lengths of 
11 mm and 2.5 mm. 
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Figure 5. Variation of inverse of compliance 1/C versus nomalized delamination length a/B for 
the ECT specimens. 
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(a) ENF Test Specimen 



(b) ENF Test Apparatus 


Figure 6. ENF test specimen and ENF test apparatus. 
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Insert 


Figure 7. Schematic of the DCB specimen. 





Figure 8. Comparison of various modes of fracture toughness for G40-800/R6376 
graphite/epoxy laminates. 
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(a) Two-Dimensional Case 



(b) Three-Dimensional Case 


Figure 1 1. VCCT for two-dimensional four-node elements and three-dimensional eight-node 
solid elements. 
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Loading Point 


X/L 


Support Point 


Figure 12. Distribution of strain energy release rate components along the delamination front 
for specimen C5 at (a/B)=0.5. 




Loading Point 


x/L 


Support Point 


Figure 13. Comparison of strain energy release rate components along the delamination front 
for different shear modulus G 23 . 
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